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Abstract. We present maps at high spatial and spectral resolution in emission lines of CCH, C-C3H2, C4H, '^CO 

I and C^**0 of the edge of the Horsehead nebula obtained with the Plateau de Bure Interferometer (PdBI). The 

. edge of the Horsehead nebula is a one-dimensional Photo-Dissociation Region (PDR) viewed almost edge-on. 

CL(' All hydrocarbons are detected at high signal-to-noise ratio in the PDR where intense emission is seen both in 

JL I the H2 ro-vibrational lines and in the PAH mid-infrared bands. C^*0 peaks farther away from the cloud edge. 

■ Our observations demonstrate that CCH, C-C3H2 and C4H are present in UV-irradiated molecular gas, with 
' abundances nearly as high as in dense, well shielded molecular cores. 

, PDR models i) need a large density gradient at the PDR edge to correctly reproduce the offset between the 

IL" ' hydrocarbons and H2 peaks and ii) fail to reproduce the hydrocarbon abundances. We propose that a new 

. 5^ ' formation path of carbon chains, in addition to gas phase chemistry, should be considered in PDRs: because of 

, intense UV-irradiation, large aromatic molecules and small carbon grains may fragment and feed the interstellar 

5_( ' medium with small carbon clusters and molecules in significant amount. 
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1. Introduction 

Due to the ISO mission, the knowledge of interstellar dust 
has significantly progressed in the past years. With its 
sensitive instruments in the mid-infrared, ISO revealed 
the spatial distribution and line profile of the Aromatic 
Infrared Bands (AIBs at 3.3, 6.2, 7.7, 8.6 and 11.3 /im 
features), which have shed light on the emission mech- 
anism and their poss ible carriers ijBoulaneer et al.(l2000t 
iRapacioli erani2004|) . However, no definite identification 
of individual species has been possible yet because the 
bands are not specific for individual molecules. The most 
likely carriers are large polycyclic aromatic hydrocar- 
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bons (PAHs) with about 50 carbon atoms ijAUain et alJ 
Il996hl: iLe Page et all l2003|) . The ubiquity of the aro- 
matic band emission in the interstellar medium has 
triggered a wealth of theoretical and laboratory work 
in the past two decades, which have led to a revi- 
sion of astrophysical models. PAHs are now suspected 
to play a major role in both the interstellar medium 
physics and chemistry. With their small size, they are 
the most efficient particles for the phot o-electric ef- 
fect fe akes & Tielens 1994; Wcingartner fc Draindl200lt 
iHabart et al.1 l200lj) . Their presence also affec t s the 
ionization balance ijFlower fc Pineau des Foretd l2003t 
'Wolfire et al. 20 03). and possibly the formation of 
H2 ( Habart'eraLll2004|l . The role of PAHs in the neutrab 
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ization of atomic ions in the diffuse interstellar medium 
has been recently re considered by iLis zt (2003), follow- 
ing previous work by *Lep p et al ,! (fl989 ). As emphasized 
soon after their discovery l|Omontlll986HLepp fc Dalgarnd 
Il988l) . PAHs also play a role in the gas chemistry: a few 
laboratory experiments, and theoretical calculations, sug- 
gest that PAHs may fragment into small carbon clusters 
and molecule s under photon impact (C?. C^k, C^H^, etc.) 
JMlin 2003; be Pa ge et all |2nn.'^lAl1a,in et a,lJ ll flflBbO: 
iLeger et al...l989.: .Scott et alJIlOOvl) . In addition, investi- 
gation of the lifetimes of interstellar PAHs implies that 
photo-dissociation may be the main li miting process for 
their life in the interstellar medium ijVerstraete et alJ 
1200 ll) . 

It is therefore appropriate to wonder whether PAHs 
could fragment continuously and feed the interstellar 
medium with small hydrocarbons and carbon clusters. 
This hypothesis is attractive for the following reasons: 



i) Cyclopropenyhdene (c-C.^H p) is widely 



in the interstellar medium jMatthews & Irvine 1985 



distributed 
1 



a) 



Matthews et al.l[l986t ICox et al...l988; .Lucas fc Liszt 
200n^ . 

Recent works have shown that the diffuse in- 
terstellar medium is more chemically active 
than previousl y thought wit h mo l ecules as 



% 

large as C3 IIGoicoeche a et al. 2004; Oka et al 

20031 lAdamkovics et ain2003^ .Roueff et al. .200 



Maier et alllioOll) and 'c-C':iH9 ijLucas fc LiszdEoO 
widely distributed. The abundance of C3 and 
C-C3H2 arc tightly connected with those of smaller 



ratios of [Cal/pg 
JCCHl/fc-C.H, ^ 



10 



40 JOka et al. 


2003") and 


jLucas & Liszl 


200a). 



in) iThorburn et al.l l|2003|) have found a correlation be- 
tween the abundance of C2 and the strength of some 
(weak) Diffuse Interstellar Bands (DIBs). 

As PAHs are present in the diffuse interstellar medium, 
could they contribute in forming both the small carbon 
clusters (C2, C3) and larger hydrocarbon molecules which 
could be the DIBs carriers? 

Unfortunately, studies of the PAHs emission bands in 
the diffuse interstellar clouds where the carbon clusters 
have been detected is extremely difficult because of the 
low column densities, and also because the bright back- 
ground star used for visible spectroscopy prohibits the 
use of sensitive IR cameras which would be saturated. 
Photo-Dissociation regions (PDRs) are the first interstel- 
lar sources in which the A IBs have beeii found and the 
PAHs hypothesis proposed l)SellgreiJll984t iLeger fc Puged 

Il984l) . It is therefore interesting to investigate the car- 
bon chemistry in these sources. iFosse et all 1I2OOOI) and 

iTevssier et al.l l|2004j) have discussed medium spatial res- 
olution (30") observations of various hydrocarbons in 
nearby PDRs. They show that CCH, C-C3H2 and C4H 
are ubiquitous in these regions, with abundances al- 
most as high as in dark, well shielded, clouds, despite 



the strong UV-radiation. iFuente et alJ l)2003j) also re- 
port high abundances of C-C3H2 in NGC 7023 and the 
Orion Bar. Heavier m olecules may be present in PDRs as 
ITevssier et all ()2004j) report a tentative detection of CgH 
in the Horsehead nebula. PDRs and diffuse clouds there- 
fore seem to share the same carbon chemistry, but because 
of their larger H2 column density and gas density, PDRs 
offer more opportunities for detecting rare species. 

ITevssier et all (|2004l) and IPuente et all l|2003l) reckon 
the presence of carbon chains is in favor of a causal link 
between small hydrocarbons and PAHs, but they lack the 
spatial resolution to conclude. In the present work, we 
present high spatial resolution observations of one source 
studied by Tevssier_fit al. (2004), the Horsehead nebula, 
obtained with the Plateau de Bure interferometer. We de- 
scribe the observations in Sec. (21 We show the interferom- 
eter maps in Sec.jSj Finally, Sec. 0] presents a comparison 
with chemical models. 



2. Observations and data reduction 

2.1. The Horsehead nebula 

The Horsehead nebula, also called Barnard 33, appears 
as a dark patch of ~ 5' extent against the bright HII 
region IC434. Emission from the gas and dust associ- 
ated with this globule has been detected froin mid- 
IR to millimet er wavelengths (^Abergel et 111120021 120031: 
ITevssier et all .2004: Pound et al. 2003). From the anal- 
ysis of the ISOCAM images, lAbergel et all ((20031) con- 
clude that the Horsehead nebula is a PDR viewed edge- 
on and illuminated by the 09.5V star crOri at a pro- 
iected distance of 0.5° (3.5 pc for a distance of 400 pc, 
lAnthonv-Twarod l(l982l) V The far-UV intensity of the in- 
cident radiation field is Go = 60 relative to the averag e 
interstellar radiation field in Draine units l|Drainelll97|i . 
The gas density, derived from the excitation of molec- 
ular lines, and from the pe netration depth of the UV- 
radiation, is a few lO'' cm~'^ l|Abergel et aLlEoO^ . From a 



combined analy 


sis of r 


iLis fc GucstciJ 


1I2OO5 


densitv.lHabart et al. 



obtain sim ilar figures for the gas 
I2004l2005ll have modeled the emis- 



sion of H2 (from narrow band images of H2 ro- vibrational 
line), PAHs and CO, and conclude that i) the gas den- 
sity follows a steep gradient at the cloud edge, rising to 
riH = 10^ cm~^ in less than 10" {i.e. 0.02 pc) and ii) this 
density gradient model is essentially a constant pressure 
model (with P = 4 x lO^Kkms^^). 

The edge of the Horsehead nebula is particularly well 
delineated by the mid-IR emission due to PAHs, with a 
bright 7.7/im-peak (hereafter named "IR peak") reaching 
25 MJy/sr at 02000 = 05''40™53.70", (S2000 = -02°28'04". 
Fig. n shows the region observed with the IRAM PdBI 
centered near the "IR peak". Two mosaics (one for hy- 
drocarbon lines and the other for the CO lines) have been 
observed. Their set-ups are detailed in Tabled 
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Phase center Number of fields 

Mosaic 1 aaooo = 05''40'"54.27' ^2000 = -02°28'00" 7 
Mosaic 2 aaooo = 05''40'"53.00= ^2000 = -02°28'00" 4 
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C-C3H2 2l,2 — lo.i 




85.339 


6.13 X 4.75 


36 


3.1x10" 


2 


Mar. 2002 ,1 


k Apr. 


2002 


C4H-I N=9-8,J= 


19/2-17/2 




85.634 


6.11 X 4.74 


36 


2.6x10" 


2 


Mar. 2002 ,1 


k Apr. 


2002 


C4H-2 N=9-8,J= 


17/2-15/2 




85.672 


6.11 X 4.74 


36 


3.4x10" 


2 


Mar. 2002 i 


k Apr. 


2002 


CCH-1 N=l-0, J 


=3/2-1/2 F= 


2-1 


87.316 


7.24 X 4.99 


54 


3.4x10" 


2 


Dec. 2002 ^ 


£ Mar. 


2003 


CCH-2 N=l-0, J 


=3/2-1/2 F= 


1-0 


87.328 


7.24 X 4.99 


54 


2.5x10" 


2 


Dec. 2002 ^ 


£ Mar. 


2003 


CCH-3 N=l-0, J 


= 1/2-1/2 F= 


1-1 


87.402 


7.24 X 4.99 


54 


3.4x10" 


2 


Dec. 2002 i 


z Mar. 


2003 


CCH-4 N=l-0, J 


= 1/2-1/2 F= 


0-1 


87.407 


7.24 X 4.99 


54 


2.3x10" 


2 


Dec. 2002 ^ 


z Mar. 


2003 


C^*0 J=2-l 






219.560 


6.54 X 4.31 


65 


9.8x10" 


2 




Mar. 


2003 


Mosaic 2 






















'2 CO J=l-0 






115.271 


5.95 X 5.00 


65 


1.2x10" 


1 




Nov. 1999 


i^CO J=2-l 






230.538 


2.97 X 2.47 


66 


1.7x10" 


1 




Nov. 1999 



° The noise values quoted here are the noises at the mosaic center (Please remember that mosaic noise is inhomogeneous due to 
primary beam correction: It steeply increases at the mosaic edges). Those noise values have been computed on 1 kms"^ velocity 
bin. 



-3°24'00" 



■3°2800" 




Fig. 1. The field of view covered when mapping small 
hydrocarbons at 3.4 mm with the Plateau de Bure 
Interferometer (PdBI) is shown as a square over this 
ESO-VLT composite image (B, V and R bands) of the 
Horsehead nebula. Each circle features the 3.4 mm pri- 
mary beam of the PdBI at one of the 7 observed positions. 
Those positions are largely oversampled at the hydrocar- 
bon wavelength (3.4 mm) to ensure simultaneous Nyquist- 
sampling at 1.4 mm used to observe C^'^O. A linear com- 
bination of the 7 pointed observation is done to obtain the 
final dirty image. 



2.2. Observations 

2.2.1. C-C3H2 and C4H 

First PdBI observations dedicated to this project were 
carried out with 6 antennas in CD configuration (base- 
line lengths from 24 to 229 m) during March-April 2002. 
The 580 MHz instantaneous IF-bandwidth allowed us to 
observe simultaneously C-C3H2 and C4II at 3 mm using 
3 different 20MIIz~wide correlator windows. One other 
window was centered on the C^^O (J=:2-l) frequency. The 
full IF bandwidth was also covered by continuum windows 
both at 3.4 and 1.4 mm. C-C3H2 and C4II were detected 
but the weather quality precluded use of 1.4 mm data. 

We observed a seven-field mosaic on a compact, hexag- 
onal pattern, with full Nyquist sampling at 1.4 mm and 
large oversampling at 3.4 mm. This mosaic, centered on 
the IR peak, was observed during about 6h on-source ob- 
serving time per configuration. The rms phase noises were 
between 15 and 40° at 3.4 mm, which introduced position 
errors < 0.5". Typical 3.4mm resolution was 6". 

2.2.2. CCH and C^^O 

As a follow-up, we carried out observations of CCH 
at PdBI with 6 antennas in CD configuration during 
December 2002 and March 2003. We used a similar cor- 
relator setup: three 20 MHz-wide windows were centered 
so as to get the four 3.4 mm hyperfine components of 
CCH; one 20 MHz-wide window was centered on the C^*0 
(J=2-l) frequency; the remaining windows were used to 
observe continuum at 3.4 and 1.4 mm. 

Exactly the same mosaic (center and field-pattern) 
and approximately the same on-source observing time per 
configuration (~ 6h) as before were used. The rms phase 
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Table 2. Calibrator fluxes in Jy. 





B0420-014 


B0607^ 


-157 


B0528+134 




3 mm 1 mm 


3 mm 


1 mm 


3 mm 1 mm 


27 11 1999 








3.5 1.4 


30.03.2002 


4.8 


2.3 






16.04.2002 


4.8 


2.5 






22.04.2002 


4.8 


2.4 






23.12.2002 


12.5 


2.6 






18.03.2003 


12.0 7.8 


2.1 


0.87 




26.03.2003 


12.8 


2.1 







noises were between 10 and 40° except during 4 hours 
in D configuration where they were between 8 and 20° 
at 3.4 mm. The data of those 4 hours have been used to 
build the C^^O map as the 1mm phase noises were then 
low enough (between 20 to 45°). We thus ended up with 
a 6" typical resolution both at 3.4 mm and 1.4 mm. Both 
CCH and C^^^O were easily mapped while no continuum 
was detected at a level of 2 mJy/beam in a 6"-bcam. 



2.2.3. ^^CO 

As part of another project (A. Abergel, private communi- 
cation), the i^CO (J=l-0) and ^^CO (J=2-l) lines were 
simultaneously observed during 6h on-source at PdBI in 
November 1999 (only 5 antennas were then available) 
in configuration C (baseline lengths from 24 to 82 m). 
The observation consisted in a 4-fields mosaic, fully sam- 
pled at 1.3 mm. The mosaic center is slightly shifted com- 
pared to the two other observations. Weather was excel- 
lent with phase noises from 3 to 5° and 6 to 10° at 2.6 mm 
and 1.3 mm, respectively. Typical resolutions were 5" at 
2.6 mm and 2.5" at 1.3 mm. 

2.2.4. Other data: H2, IS0-LW2 and 1.2mm dust 
continuum 

The H2 v=l-0 S(l) map shown here is a small part 
of Horsehead observations obtained at the NTT us- 
ing SOFI. The resolution is ~ 1". Extensive explana- 
tions o f data reduction an d analysis are discussed else- 
where llHabart et all I2004L 12005). The IS0-LW2 map 

Ilk 



(pubhshed by I Abergel et al]T2003D ) shows aromatic fea- 
tures at 7.7^m with a resolution of ~ 6". The 1.2 mm dust 
continuum has been obtained at the IRAM-30m telescope 
wi th a resolution of ~ 1 1" and has already been presented 
bv lTevssier et al.l I 



2.3. PdBI data processing 

All data reduction was done with the GILDAS^ softwares 
supported at IRAM. Standard calibration methods using 
close calibrators were applied to all the PdBI data. The 



calibrator fiuxes used for the absolute flux calibration are 
summarized in Table El 

Following ICueth et ahl lll996l) . single-dish, fully 
sample d maps o btained at the IR AM-30m tele- 
scope (|Tevssier et al. 20o| lAbergel et al. 2003)) were 
used to produce the short-spacing visibilities flltered out 
by every mm- interferometer (e.g. spatial frequencies be- 
tween and 15 m for PdBI). Those pseudo- visibilities were 
merged with the observed, interferometric ones. Each mo- 
saic fleld were then imaged and a dirty mosaic was built 
combining those fields in the following op timal way in 
terms of signal-to-noise ratio llrxuethll200lfl : 



E, 



^ See http://www.irain.fr/IRAMFR/GILDAS for more infor- 
mation about the GILDAS softwares. 



In this equation, J{a,S) is the brightness distribution in 
the dirty mosaic image, the Bi are the response functions 
of the primary antenna beams, the Fi are the brightness 
distributions of the individual dirty maps, and the Ci are 
the corresponding noise values. As may be seen on this 
equation, the dirty intensity distribution is corrected for 
primary beam attenuation. This implies that noise is inho- 
mogeneous. In particular, noise strongly increases near the 
edges of the field of view. To limit this effect, both the pri- 
mary beams used in the above formula and the resulting 
dirty mosaics are truncated. Standard level of truncation 
is set to 20% of the maximum in GILDAS. In our case, 
the intensity distribution does not drop to zero at all field 
edges. Hence, we used a much lower level of truncation 
of the beam {i.e. 5%) to ensure a better deconvolution of 
the side lobes of the sources sitting just at the field edges. 
We then use the standard adaptation to mosaics of th e 
Hogbom CLEAN algorithm to deconvolve l)Guethl 120011) . 
The sharp edge of the II2 emission defines a boundary 
that may be used as an a priori knowledge in the decon- 
volution of the PdBI images: We use this boundary as a 
numerical support (in the language of signal processing) 
to forbid the search of CLEAN components outside the 
PDR front {i.e. in the direction of the exciting star). We 
finally truncated the noisy clean mosaic edges using the 
standard truncation level. It may be noted that the C4H 
maps are particularly difficult to deconvolve due to their 
low signal-to-noise ratio, S/N < 10 to 15. 

3. Results 

3.1. Maps 

The PdBI maps are shown in Fig. [21 and P tog ether 
with the 7^m ISOCAM im age llAbergel et allEoO^ . the 
1.2 mm dust emission map (•Tevssier et al."2004') and the 
map of the II2 2.1 /im line emission (,Habart et al. 2004L 
I2OO5I) for comparison. For all lines, we obtained excel- 
lent spatial resolutions, similar to or even better than the 
ISOCAM pixel size of 6" (see Table P). Fig. El shows the 
maps in the natural Equatorial coordinate system while 
Fig.Oshows the maps in a coordinate system where the x- 
axis is in the direction of the exciting star and the y-axis 
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Fig. 2. InteKrated emission maps obtained with the Platea u de B ure Interferometer. Maps of i) the H2 v=l-0 S(l) 
emission (Haba rt et al. 2004. 20Q51. ii ) the mid-IR emission l)Abergel et al..2Q 03. labeled IS0-LW2) and Hi) the 1.2 mm 
dust contimuim~ ljTe"vssier et al.ll2004l labeled 1.2mm) are also shown for comparison. The center of all maps has been 
set to the mosaic 1 phase center: RA(2000) = 05h40m54.27s, Dec(2000) = -02°28'00". The map size is 110" x 110", 
with ticks drawn every 10". Either the synthesized beam or the single dish beam is plotted in the bottom, left corner. 
The emission of all the lines observed at PdBI is integrated between 10.1 and 11.1 kms^^. Values of contour level are 
shown on each image wedge (contours of the H2 image have been computed on an image smoothed to 5" resolution). 
The sharp edge of the H2 emission (upper right panel) defines a boundary, which is used as a numerical support (in 
the language of signal processing) for deconvolution of the other images. This deconvolution support is overplotted in 
red on every panel. 



defines an empirical FDR edge that corresponds to the 
sharp boundary of the H2 emission (i.e. the maps have 
been rotated by 14° counter-clockwise and horizontally 
shifted by 20"). The latter presentation enables a much 
better comparison of the FDR stratification. 



The main structure in all hydrocarbon maps is an 
approximatively N-S filament, following nicely the cloud 
edge and corresponding closely to the mid-IR filament on 
the IS0-LW2 image. A weaker and more extended emis- 
sion is also detected, which has no counterpart in the 
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Fig. 3. Same as Fig.|21except that maps have been rotated of 14° counter-clockwise around the image center to bring 
the exciting star direction in the horizontal direction as this eases the comparison of the PDR tracer stratifications. 
Maps have also been horizontally shifted by 20" to set the horizontal zero at the PDR edge delineated as the vertical 
red line. Horizontal red lines delimit the two lanes that have been vertically averaged to produce the two series of cuts 
shown in Fig. El 



IS0-LW2 image and can be attributed to the bulk of the 
cloud. It is interesting to note that the hydrocarbon emis- 
sion presents a minimum behind the main filament, and 
a weaker secondary maximum within the extended emis- 
sion. The hydrocarbon emission is stronger on the edges 
of the dust 1.2 mm emission and avoids the region of max- 
imum dust emission where the gas is likely denser. This 
confirms a tendency revea led bv chemical sury evs of dense 
co res (study of TMC-1 b ylPratar) et alJ l|l997l) and L134N 
bv lDickens et alJ l)200(]|) :' lFossel l)2003l) ): i.e. carbon chains 



(CCH, C4H,...) generally avoid the densest and more de- 
pleted cores. 

Even at the high spatial resolution provided by the 
plateau de Bure Interferometer, the maps of all hydro- 
carbons remain very similar. Detailed inspection of the 
maps shows small differences between CCH and C-C3H2, 
but these do not affect the overall similarity. Indeed, the 
joint histogram describing the correlation of line maps 
for i) the two most intense CCH lines, ii) C-C3H2 and 
CCH, and iii) C4H and CCH are displayed in Fig. ^ 
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-15" < 5y < 0" ^ 0" < 5y < 15" 




50 -50 50 -50 

6x (arcsec) (5x (arcsec) 



Fig. 5. Emission profiles along the exciting star direction (PA — -104° in equatorial coordinate system). To improve 
signal-to-noise ratio, those emission profiles have been integrated along the perpendicular direction between —15" < 
6y < 0" (left column) and 0" < Sy < +15" (right column). The comparison of those two series of mean cuts gives an 
idea of the influence of the local conditions (either excitation effects or chemical differentiation) . We show from top to 
bottom, H2 v=l-0 S(l) (full resolution in black, smoothed at a 5"-resolution in red), ^^CO J=2-l, ^^CO J=l-0, C^^O 
3—2-1, CCH lines (1 black, 2 green, 3 red and 4 blue, respectively following the ratios 1:0.5:0.5:0.25), C-C3H2, C4H 
lines (1 black, 2 red), IS0-LW2 and the 1.2mm dust continuum. For the PdBI data, the 3-cr noise level is indicated 
by the dashed lines. It rises at the cut edges due to the primary beam correction. Note that the fields of view of the 
^^CO and C^^O data are smaller than the field of view of the hydrocarbon data because of the smaller mosaic size 
and/or the higher frequency. Finally, the blue vertical lines give guideline to localize the hydrocarbon peaks. They 
have been drawn at Sx = +12" (left column) and Sx = +15" (right column). Note that the bumps at Sx = —20" in 
the IS0-LW2 cut (left column) and at Sx — 70" in the H2 cut (right column) are due to field stars. 
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Fig. 4. Joint histogram of the integrated emission oi i) 
the second brightest CCH line (top), ii) C-C3H2 (middle) 
and Hi) one C4H line (bottom) vs. the main CCH line. 
The value at a given position of this joint histogram is the 
percentage of pixels of the input images whose intensities 
lies in the respective vertical and horizontal bins. Only im- 
age pixels lying inside the deconvolution support (shown 
in Fig. 121 have been used in the histogram computation. 
Contour levels are set to 0.125, 0.25, 0.5, 1, 2, 4 and 8% 
of points per pixel. 



As expected the two CCH lines are extremely well corre- 
lated as illustrated by the elongated shape (approaching 
a straight line) of the joint histogram. The correlations 
between C-C3H2 and CCH, and between C4H and CCH 
are excellent too, though the signal-to-noise ratio is not 



as good for C4H. For this plot, we have used all points 
lying inside the support used for the deconvolution. 

The high resolution C-C3H2 map appears to show more 
structure than the CCH maps, particularly in the well 
shielded cloud interior (on the left hand side of the main 
filament). This effect seems real since it does not appear 
for the satellite CCH line maps, which have similar inten- 
sities and signal-to-noise ratio as the C-C3H2 map. The 
C4H maps are too noisy for a detailed analysis but are 
nevertheless very well correlated with the CCH map. The 
corre lations found at low spatial resolution (Teussi^r et al\ 
\p.nnA) are not an artifact hut persist at high spatial reso- 
lution. 

The correspondence of hydrocarbons with CO and 
C^^O is not as good. The C^^O (J=2-l) map presents two 
maxima, located on either side of the CCH peak along a 
N-S direction: The CCH peak is associated with a local 
minimum of C^^O emission. Also, the C^^O emission peak 
is displaced farther inside the cloud (East) compared to 
CCH and the other hydrocarbons. 

To illustrate further the differences in the spatial dis- 
tribution of CO, C^^O and the hydrocarbons, we show 
two series of cuts across the PDR in Fig. El The UV- 
radiation comes from a Ori far to the right side of the 
Fig. 13 The cuts have been taken along the a Ori direction 
[i.e. PA= —104°). The main peak for all hydrocarbons is 
located near an offset of Jcc ~ 12 — 15" at less than 5" 
of the H2 peak. The IS0-LW2 peak is located half-way 
between hydrocarbons and H2 peaks. Intense '^^CO emis- 
sion, in both the J=l-0 and J=2-l lines is also detected 
in the same region, while the C^^O (J=2-l) emission rises 
farther (at least 5") inside the cloud. 

As shown in Fig.|Hl the ^^CO (J=2-l) emission (con- 
volved at the same angular resolution as the ^^CO J=l-0 
transition) is very bright (> 50 K at 10.6 km s^^, the line 
peak velocity) and more intense than ^^CO (J=l-0) in the 
most external layers of the PDRs, facing directly crOri. 
The line intensity ratio Tb(l — 0)/r;,(2 — 1) rises from 
~ 0.3 to ^ 0.8 from West to East. Combined with the high 
brightness temperature detected for both lines, the higher 
brightness temperature of the ^^C0(2-l) line is a clear sign 
of the presence of warm and dense gas. We have estimated 
the kinetic temperature using an LVG model. We assumed 
that the emission is resolved and fills the beam. We ex- 
plored the kinetic temperature dependence upon the den- 
sity by solving for 5 different proton densities going from 
1.6 X lO-^cm-^ to 10^ cm . Under these hypotheses, the 
^^CO line intensity ratio and brightness temperature con- 
strain the kinetic temperature to increase from 60 K in the 
inner PDR (15" = 0.03 pc from the PDR edge) to more 
than 100 K in the outer layers for proton densities larger 
than 4 x lO^cni^^. For lower proton densities, the kinetic 
temperature still starts from 60 K in the inner PDR but 
increases much more stiffly. The k inetic temperature d e- 
rived from single dish observations ijAbergel et alJ2003|) is 
lower, in the 30—40 K range and corresponds to the bulk of 
the cloud, rather than to the warm UV-illuminated edge. 
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Fig. 7. CO spectra (convolved at the same angular resolution) along the direction of the exciting star at dy — —2.5". 
In those cuts, the label Sx — 13" indicates the IR peak position (cf. table [SJ. Note that ^■^CO J=2-l peak intensity 
decreases at positions Sx = 13" and 15" while "'^^CO J=l-0 increases even reaching its maximum at Sx — 15". In 
addition, both ^^CO lines show a small but clear dip {i.e. the center channel intensity is lower than its first neighbours) 
at Sx = 15". Finally, while C-'^^O spectra are very close to Gaussians, ^^CO spectra show asymmetric profiles. Spectra 
cuts at other close Sy values show the same trends. 



3.2. Abundances 

We have computed the CO and hydrocarbons column 
densities at three representative positions in the maps: 
the "IR peak" where the PAH and hydrocarbon emis- 
sion is the largest, the "IR edge" 10" West which repre- 
sents the region with the most intense UV-radiation and 
a "Cloud" position behind the IR filament. Table 13 hsts 
the derived column densities and abundances relative to 
the total number of protons for these 3 positions. We have 
used a LVG model with different uniform total hydrogen 
density^ (from 10'*cm~'^ to 10^ cm^'^) and a kinetic tem- 
perature of 40 K for the "cloud" position, and between 60 
and 100 K for the IR positions. The variance of the col- 
umn densities therefore reflects both the systematic effect 
due to the imperfect knowledge of the physical conditions, 
and the random noise of the data. In most cases, the for- 
mer contribution is the largest. The H2 column densities 
are derived from the dust 1.2 mm emission assuming the 
same dust properties for all positions but a dust temper- 
ature range of 20 to 40 K for the "Cloud" position and 40 
to 80 K for the IR positions. 

The LVG solution implies a typical ^^CO column den- 
sity of 2 X 10^^ cm^^. This is inconsistent with the derived 
column density of C^^O and the local IS M ^^O/^^O el- 
emental ratio ('560. IWilson fc RoodlHool . Fig. □ shows 



"Total hydrogen density" is a short cut for the total density 
of hydrogen in all forms. 



clear indications of self-absorption of the ^^CO spectra 
(asymmetries and dips in the top of the line profiles) while 
the C^^O spectra are Gaussian. The same behaviou r is 
seen i n the single dish data discussed by Abcrgcl et alJ 
l|2003l) (cf. their Fig. 5). This explains why the LVG solu- 
tion does not succeed in correctly inferring the ^^CO col- 
umn density. Conversely, the C^'^O abundance relative to 
H is fairly constant for all positions at [C^^O] = 1.0 x 10^^. 
Assuming a local ISM -"^^0/^*0 elemental ratio, this cor- 
responds to a CO abundance relative to the total num- 
ber of hydrogen atoms of [CO] — 5.6 x 10^^, in rather 
good agreement with the gas phase abundance of car- 
bon derived from CO in warm molecul ar clouds, and to 
the carbon abundan ce in diffuse clouds ijLacv et al.lll994 
ISofia fc Me'veill200l|) . In addition, using IRAM-30m spec- 
tra of "CO and C^^O published by Abergel et al. (20q3), 
we found ["CO]/[C^^O] 7. This good agreement with 
the local ISM isotopic ratio make us confident that we can 
use our LVG analysis on the PdB I C^^O spectra to esti - 
mate the CO density. According to lLis fc GuestenI l)2005|) . 
atomic carbon is less abundant than CO in the PDR. The 
peak column density of neutral carbon, observed with a 
15" beam, is ^ 1.6 x 10^^ cm^^ corresponding to a carbon 
abundance of [C] = 5 x 10^^. Even if we take into ac- 
count the difference in linear resolution, we don't expect 
an increase of the column density larger than a factor of 
two based on the comparison of the low resolution single 
dish data with the interferometer maps of other tracers. 
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Fig. 6. Spatial variation along the direction of the exciting 
star of a) the ^^CO J=2~l brightness temperature (con- 
volved at the same angular resolution as the J=l-0 transi- 
tion), b) 1-0/2-1 ratio and c,d,e) the kinetic temperature. 
{6x,5y) offsets refer to the coordinate system defined in 
Fig. 13 The cuts at dy = +10" are drawn in blue, those at 
6y = +5" in green and those Sy = +0" in red. The data 
cuts have been taken from the 10.6 km s^^ velocity channel 
corresponding to the ^^CO line peak. The kinetic temper- 
ature is derived from an LVG model assuming i) unity 
beam filling factor and ii) uniform total hydrogen density 
of nn = 1-6 X 10"* cm"-^ (dashed lines), 2 x 10^ cm"^ (full 
lines), 4 x lO^'cm^'^ (dotted lines) and 10^ cm^^ (dotted- 
dashed lines). Arrows indicate lower limits. 



Table 3. Molecular column densities and abundances at 3 
different positions of the PDR named "Cloud" , "IR peak" 
and "IR edge". Equatorial offsets refer to the Mosaic 2 
map center given in Table Q (^2:, Sy) offsets refer to the 
coordinate system defined in Fig. |3 H2 column densities 
have been derived from the 1.2 mm dust continuum emis- 
sion using a dust temperature range of 20 to 40 K for the 
"Cloud" position and 40 to 80 K for the IR positions. 
Others column densities used LVG models with a rep- 
resentative set of densities and kinetic temperature. 1-cr 
uncertainties thus reflect the systematics due to the ap- 
proximate knowledge of density and kinetic temperature. 
Abundances are computed with respect to the number of 
protons, i.e. [X] = 0.5 A^(X)/iV(H2). 





SRA 


SDec 


5x 


Sy 


Cloud 


+6" 


-4" 


+24.9" 


-5.3" 


IR peak 


-6" 


-4" 


+ 12.2" 


-2.4" 


IR edge 


-12" 


-4" 


+7.4" 


-1.0" 





Quantity 


Unit 




Cloud 


IR peak 


IR edge 




Si. 2 mm 


mJy/Beam 


38+2 


35+2 


12+2 




iV(H2) 


10^^ cm" 


2 


27+9.5 


10.5+4 


3.6+1.7 






10^^ cm" 


2 


5.8+0.5 


4+0.5 


1+0.3 




iV(CCH) 


10" cm" 


2 


5.5+1 


30+5 


11+3 




iV(c-C3H2) 


10^^ cm" 


2 


2.3+0.7 


24+10 


9.5+5 




N{C4'R) 


10^^ cm- 


2 


20+10 


40+10 


37+10 








Quantity 


Unit 




Cloud IR peak IR edge 






10-'^ 




1.07 


1.9 


1.4 




[CCH] 


10-** 




0.10 


1.4 


1.5 




[C-C3H2] 






0.43 


11.4 


13.2 




[C4ff 


10-3 




0.37 


1.9 


5.2 



Finally, though the H2 column densities are fairly similar 
at the "IR peak" and "cloud" positions, the abundances 
of hydrocarbons are larger by a factor of at least 5.0 at 
the "IR peak" . The abundances seem to be even larger at 
the "IR edge" than at the cloud position. 

4. Discussion 

4.1. Comparison with models 

We have used a mo nodimen sional 

PDR code ijLe Petit et alJ I2OO2L 

http://aristote.obspm.fr/MIS/) for modeling the 
observations of the Horsehead nebula. The slab geometry 
is locally appropriate as seen on Fig. O We did not take 
into account projection effects as the sourc e is viewed 
almost edge-on. Indeed, iHabart et alJ l)2005(l show that 
the main effect of the PDR possible small inclination 
(< 6°) is to enlarge the peak profiles and to shift them all 
compared to the model edge. The model includes a de- 
tailed treatment of the photo-dissociation of H2 and the 
CO isotopes as well as the statistical equilibrium of their 
rovibrational (rotational, respectively) states in a steady 
state approach. The parameters of the model include the 
elemental abundances, the cosmic ray ionization rate, 
the scaling factor of the interstellar ultraviolet radiation 
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Fig. 8. Predictions of the spatial variation of the abundance relative to H2, using a unidimensional FDR code. For 
each model, the abundance of the population of the the upper level of the 2.12 /xm H2 line {i.e. v=l, J=3), written 
[H2], is shown in linear scale (top). The C, CO and hydrocarbon abundances are shown in logarithmic scale (bottom). 
The modeled cloud is illuminated from the right-hand side. The Jx-axis origin has been set so that [Hj] peaks at the 
position of the observed H2 peak {i.e. Sx = 10"). The vertical dotted blue hue indicates the peak of the C-C3H2 and 
C4H abundances. Each model is described in 3 hues: i) the density structure, ii) the UV-field properties and Hi) the 
chemical network used. 6 different models are compared here. Our reference model is labeled A. The total hydrogen 
density is kept uniform at a value of tt-h = 10^ cm~^. The far-UV intensity of the radiation field is Go = 100 (in Draine 
units) and the extinction curve is the mean Galactic one. The chemical network rate file is the New Standard Model 
one with minor modifications described in the text. Charge exchange reactions between C+ and FAHs are not taken 
into account. The gaseous sulfur abundance is low compared to solar {i.e. S/H = 5.8 x 10~*). The parameters varied 
in other models are emphasized in red. Model B uses a different extinction curve. Model C adds reactions of charge 
exchange between C+ and FAHs. Model D decreases the un iform total hydrogen den sity. Models E and F use the 
density profile derived from the model of the H2 observations IjHa 
sulfur abundance {i.e. S/H = 10~^). 



h ydrogen den s 



Model F uses a solar gaseous 
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field (ISRF) measured in Draine units, the density profile 
and the grain parameters. As the observations involve 
complex carbon molecules, we have used the so called 
"new standar d model" chem ical rate file of Herbst and 
collaborators llLee et al ]|1998^. available on the web site^ . 
In a previous paper l|Tevssier et "ani2004j) . we have found 
that the other ext ensive chemical rate file provided by 
the UMIST group l)Le Teuff et al.ll2000() gave close results 
on the carbon chain molecules. As C'^^O observations 
are reported, we have added to this reaction set, the 
main isotopic molecules involving a nd introduced 
the c orresponding fractionation reactions ijGraedel et alJ 
We have also introduc ed the photo-dissociation 
rates given by Ivan Dishoeckl (|^988), when available, 
which have been calculated specifically with the Draine 
ISRF and which were different from the values reported 
in the chemical rate file. The resulting chemical network 
involves about 450 chemical species and 5000 reactions. 
Only the most stable isomeric forms of hydrocarbons are 
considered here. 

We define a reference model (hereafter named model 
A) of the Horsehead nebula as a uniform sheet of gas and 
dust of total hydrogen density rin = 10^ cm"'^ exposed 
to a ISRF of 100 measured in Draine units. The cosmic 
ray ionization rate has a value of 5 x 10^^^ s^^ and the 
elemental abundances are as follows: C/H = 1.38 x 10^"*, 
0/H = 3.02 X 10-4, iSQ/H = 6 x lO"'^, N/H = 7.95 x 10-^ 
S/H = 5.8 X 10-^ Cl/H = 1.86 x 10"^, P/H = 9.3 x 10'^°, 
Fe/H = 1.7 X 10-^ Mg/H = 10"^, Na/H = 2.3 x 10"^. 
T he properties o f the g rains are the same as described 
in iLe Petit et alj l| 2002t). i.e. t he size distribution law is 
taken from lMathiset^ with an exponent of —3.5 

and we describe the attenuation of grains from the far- 
UV to the visible via the galactic extinction curve given 
as an anal ytic function of 1 / A by includ ing: the coefficients 
derived bv iFitzpatrick fc Massal 1 1988() . Charge exchange 
reactions between 0+ and PAHs are not taken into ac- 
count. The gas to dust mass ratio is 100. 

Fig. IHl shows i) the abundance of the H2 rovibra- 
tionnally excited in the v=l, J=3 level at the origin of 
the 2.12 /zm line (this abundance is hereafter referred to 
as [H2]) and ii) the C, CO and hydrocarbon abundances 
for this reference model and 5 variants. We ensure that 
the [H2] peak position is set at Sx — 10" as in the ob- 
servations. Our reference model correctly reproduces the 
observed 3 to 5" offset between the hydrocarbon and II2 
peaks. The C^^O also peaks behind the hydrocarbons at 
Sx = 20 — 25". However, the H2 profile is not correctly 
modeled here. 

In model B, we replaced the Galactic extinction 
curve by one more representative of molecular gas. 
We have chosen HD 147889 in Ophiuchus. Its extinc- 
tion curve has a rather s trong far-UV rise {Eb-v — 
1.09 IFitzpatrick fc Mas'sal (|l98i)). Its ratio between the 
total and selective extinctions, Ry, is 4.2 a typical fig- 
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Fig. 9. Spatial variation of the total hydrogen density in 
models A to F. In model E and F, the density is increasing 
as a power law of scaling exponent 4 on the first 10" and 
then is kept constant at a value of 2 x 10^ cm"^. As in 
Fig. IHl the X-axis origin has been set so that [Hj] peaks 
at the position of the observed 2.12 /im, H2 line peak (i.e. 
Sx = 10"). 



http : //www. physics . ohio-state . edu/ "eric/research, 
files/ cddata. julyOS 



ure for molecular gas (jCordon et aljlioost ICardelli et alJ 
1989). The PDR stratification does not qualitatively 
change compared to model A: It is just compressed. In 
model C, we added reactions of charge exchange between 
C"*" and PAHs. This enhances the neutral atomic carbon 
abundance but does not have a large effect on the hydro- 
carbons: Only CCH peaks closer from the H2 peak com- 
pared to model A. Neither model B nor C improves the 
modeling of the H2 profile. 

As shown by model D, E and F, the density structure 
has a major impact on the PDR structure. Fig. |5| shows 
the density profiles associated to each model. When keep- 
ing the total hydrogen density uniform but decreasing its 
value to 2 X 10* cm"^ (as in model D), the carbon and hy- 
drocarbon abundance peaks highly broadened and shifted 
inward by more than 20", a prediction clearly violated by 
the high resolution PdB I data. Model E and F uses a den- 
sity profile provided by iHabart et aD l|2004l l2005l) to fit 
the 2.12(Um-H2 emission. Indeed, the [Hj] profile qualita- 
tively changes (it is now a peak rising from zero at the 
PDR edge) but it also reproduces the H2 filament width. 
Those two models which impose a steep total hydrogen 
density gradient at the PDR edge, are the only ones that 
succeed to correctly reproduce the offset between the hy- 
drocarbon and H2 peaks as well as the form of the H2 
peak. The sole difference between models E and F is the 
gaseous sulfur abundance: sulfur is depleted from the gas 
phase in model E (S/H = 5.8 x 10"^) while the gaseous 
sulfur abundance is solar in model F (S/H = 10~^). 
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Fig. 10. Comparison between our two best models 
(curves) and observed (points with error bars) abundances 
of tlie small hydrocarbons. CCH is shown as a green, plain 
line, C-C3H2 as a dashed, red line and C4H as a blue dotted 
line. The top and bottom panels respectively show abun- 
dances relative to the total hydrogen density and CCH. 
The dashed vertical line shows the position where the to- 
tal hydrogen density profile swaps from a steep gradient 
to a constant. 

Fig. EH is a zoom for our two best models {i.e. E and 
F) of the spatial variations of the abundances of hydro- 
carbons relative to i) total hydrogen density (top panel) 
and ii) CCH (bottom panel). The observed abundances 
are overplotted with their error bars. The dashed verti- 
cal line separates the zone where the proton gas density 
is constant from the zone where the proton gas density 



rapidly decreases outward. This latter zone is associated 
with the PDR. The sulfur elemental abundance has differ- 
ent effects in those two regions. In the region of moderate 
visual extinction {i.e. the"IR edge" and the "IR peak" 
where Ay ^1), the charge transfer reaction between C+ 
and S leading to S"*" and C reinforces the abundance of 
neutral carbon and thus enables the formation of car- 
bon chains via the rapid neutral-carbon atom reactions. 
However this effect is small. Indeed this is in the dark 
region where the sulfur elemental abundance has a large 
effect. When the sulfur abundance is solar, the small car- 
bon chains C2, CCH, C2H2, C3H, C3H2 and C4H react 
with S+ to give C2S+, CCS+, HC2S+, C3S+, HC3S+ and 
C4S+. In this main destruction path of the small carbon 
chains, one hydrogen atom is released impairing the ref- 
ormation of the carbon chains. When S is higly depleted 
as in Model E, this destruction mechanisms is superseded 
by other paths involving C^. Those paths form carbon 
chain ions which in turn contributes to the formation of 
other carbon chains. Overall, model E {i.e. low S/H) per- 
forms better in the comparison with observed abundances. 
The only exception is the n(c-C3H2)/n(CCH) ratio at the 
"cloud" position. We will thus use model E only for com- 
parison with the observations. At the IR peak (median 
point at Sx — 15.5"), CCH abundance is correctly repro- 
duced while C-C3H2 and C4H abundances are underes- 
timated by at least a factor of 3. Discrepancies are much 
higher both at the "cloud" (point to the right at Sx = 27") 
and the "IR edge" (point to the left at 6x = 9") posi- 
tions. In the UV-illuminated edge, the modeled [CCH] 
has a quite shallow increase with Sx while the modeled 
[C-C3H2] and [C4H] share the same steep abundance pro- 
file. By contrast, the observed ("IR edge") abundances 
are very similar for the 3 species, reflecting the very good 
spatial correlation between the different hydrocarbons (see 
Fig. 0)) . This discrepancy is independent of our knowledge 
of the total hydrogen density as it is also seen when com- 
paring abundances relative to CCH. 

In summary, none of our models is able to correctly 
reproduce the relative stratification of H2 and small hy- 
drocarbons. Comparison of model A and C shows that 
to reproduce the observed offset between hydrocarbon 
and H2 peaks, we need a high total hydrogen density 
(10^ cm^^). By varying the profile density (model E and 
F), a shallow total hydrogen density increase at the PDR 
edge is needed for reproducing the profile of the 2.12 /im 
H2 line. However, the shallower the total hydrogen density 
increase, the larger the modeled offset between H2 and the 
hydrocarbons. A good compromise is provided by a total 
hydrogen density profile increasing as a power law with a 
scaling exponent 4 on the first 10" a nd then constant at a 
value of 2 X 10^ cm~^. lHabart et al"l| 20Q5) shows that this 
model essentially correspond to a constant pressure model 
(with P = 4 X lO^Kkms"^). Finally, model E with low 
sulfur elemental abundance performs better than model F 
with solar abundance. Nonetheless, even model E do not 
succeed to reproduce the good hydrocarbon correlation 
seen in the illuminated part of the PDR: while CCH is 
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correctly predicted to have a smooth abundance increase, 
modeled C-C3H2 and C4H abundances show a much too 
steep increase. 

4.2. Can the fragmentation of PAHs contribute to the 
synthesis of small hydrocarbons? 

Examining the model predictions in more details, three 
hypotheses can be proposed for explaining the discrepan- 
cies between model calculations and observations: 

i) The photo-dissociation rates used in the models may 
be incorrect. As the main destruction process near 
the cloud edge is photo-dissociation, the actual val- 
ues of the photo-dissociation rates are critical for an 
accurate prediction. However, similar results are ob- 
tained with the UMIST95 and NSM rate files. The 
photo-dissociation rates for C-C3H2, C3H and C4H are 
10~^ for both rate files, and differ by a factor 
of two for CCH {i.e. 0.51 x 10"^ s^^ for UMIST95 
and 10~^ s~^ for NSM). The photo-dissociation rates 
of larger chains are similar. In most cases, except for 
CCH and acetylene, the numbers given i n the rate files 
are n ot well documented. For instance, Ivan Dishoeckl 
l)l988il discusses the photo-dissociation rate of C-C3H2 
and concludes that it is accurate within an order of 
magnitude. More accurate photo-dissociation rates are 
clearly needed for the carbon chains and cycles. Recent 
calculations have been performed for C4H showing 
that the photo-dissociation threshold is 5.74 eV, but 
that efficient photo-dissociation requires more en- 
erget ic photons, typically above 6.5 eV l|Craf et al.l 
l2onih . However, it is unlikely that the rates are low 
enough to explain the large discrepancies between the 
models and the data since these molecules are known 
to be sensitive t o UV-radiation I Jackson et al.lll99ll 
ISnrig et ^11199411 . 
ii) Another possibility is that the chemical networks are 
missing important reactions for the synthesis of hydro- 
carbons. Neutral-neutral reactions are progressively 
included in the rate files, but still are much less nu- 
merous than ion-molecule reactions. It is now known 
that atomic carbon, dia tomic carbon and CCH may re- 
act w i th hydrocarbons llKaiser et al ] l20n.1:IStahl et alJ 
l2002t iMebel fc Raised l2002|) . More work remains to 
be done. However, preliminary tests using a more ex- 
tended data base of chemical reactions have not led to 
significant improvement. 
Hi) The excellent spatial correlation between the mid-IR 
emission due to PAHs and the distribution of car- 
bon chains suggests a last hypothesis: the fragmen- 
tation of PAHs due to the intense far UV-radiation 
could seed the interstellar medium with a varietv of 
carbon clusters, chains and rings fScott et al.' "1997', 
Verstraete etall l200lt iLe Page et al. .2003; .Joblir 



2003i and references therein) . These species would then 
further react with gas phase species (C, C+, H, H2, 
etc.) and participate in the synthesis of the observed 



hydrocarbons. iFuente et al.1 l)2003j) also favor this ex- 
planation for explaining the abundance of C-C3H2 they 
observed in other PDRs. 

A correct exploration of this third hypothesis needs a 
good qualitative and quantitative description of both the 
fragmentation and reformation of PAHs, which is out of 
the s cope of t his p aper. We here give only a few indica- 
tions. |l986) pioneered the attempts to understand 
the role of P AHs in the interstellar ch emistry Elaborating 
on this work. lLepp fc Dalgarnol l)l988|) suggested that par- 
ticipations of PAHs in the ion chemistry of dense clouds 
lead to large increase in the abundances of small hydro- 
carbons. Indeed, when the PAH fractional abundance ex- 
ceeds ~ 10~^, the formation of PAH^ triggers mutual 
neutralization of the positive atomic and molecular ions 
and introduces new pathways for the formation of complex 
molecules. The equilibrium abundances of neutral atomic 
carbon C, CCH and C-C3H2 may thus be enhanced by two 
order of magnitude. By comparison, our model C which 
includes charge exchange between C+ and PAHs, shows 
a decrease of C4H and C-C3H2 abundances by at least 
an order of magnitude in the dark region. Introduction 
of mutual neutralization between C+ and PAH^ could be 
an interesting alternative to our "artificial" lowering of the 
sulfur abundance. We are currently acquiring CS data at 
PdBI to constrain the S chemistry independently. 

iLepp et all l|l988fl suggested that the ion chemistry 
of diffuse clouds has little impact on the CH, OH and 
HD abundance, but can lead to large increase in the 
abundance of other species (H2, NH3 and most notice- 
ably CH4 and C2H2) by successive reactio ns of PAH 
and P AH~ with carbon and hydrogen atoms. iTalbi et alJ 
suggested that Coulombic explosion of doubly ion- 
ized PAH could create C-C3H2 through the electronic dis- 
soc iative recombin ation of C3Hj^. Laboratory experiments 
by I.Tochims et alJ l|1994) suggested that PAHs with less 
than 30-40 carbon atoms will be UV-photodissociated in 
Hi regions while lar ger ones will be stable. Based on those 
results, models by lAllain et al.l l)l996blafl indicate that 
only PAHs with more than 50 carbon atoms survive the 
high UV radiation field of the diffuse interstellar medium, 
whereas smaller PAHs such as coronene or ovalene are 
destroyed by the loss of acetylenic groups. Destruction 
timescales are a few years for neutral species and typi- 
cally five time shorter for the corresponding cations. All 
those reactions starts from neutral or cation PAHs. They 
will be in competition with charge exchange and mutual 
neutralization discussed above. Mutual neutralization has 
a maximal effect in the transition region where the gas is 
molecular but the electronic abundance is significant. This 
region corresponds more or less to the region of maximum 
emission from the PAHs or slightly deeper in the molec- 
ular cloud. All other cited reactions are more efficient to- 
ward the illuminated edge whe re PAHs are mainly neutral. 
Recently. iLe Page et al.l l|2003|) discussed the possibility of 
addition reactions with ionized carbon, starting from the 
high reaction rate between C+ and anthracene measured 
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by ICanosa et alJ 1|1995|) . If similar reaction rates persist 
for heavier PAHs, addition reaction with carbon would be 
very efficient for counteracting the destruction by far-UV 
photons. 

From the observational point of view, the mid-IR emis- 
sion due to PAHs is extended in interstellar clouds. On 
the other hand, a detailed analy sis of the mid and far- 
IR i mages obtained bv IR AS led iBoulanger et al. (199(1) 
and iBernard et al.l l|l99,'^ to conclude that PAHs disap- 
pear in the dense cold cloud int eriors, probab l y bec ause 
they coagulate and/or condense. IStepnik et al.l l)2003|) de- 
scribe a convincing case for such a process in a sm all fil- 
ament of the Taurus cloud. iRapacioli et alJ l)2004) have 
found clear evidences for spatial variations of the aromatic 
infrared band profiles, likely due to the spatial variation 
of the nature of their carriers. A sophisticated analysis of 
ISOCAM-CVF data allow them to separate the mid-IR 
spectra of the ionized and neutral PAHs from the spec- 
tra of carbonaceous very small grains (possibly PAHs ag- 
gregates). The very small grains are located at a larger 
distance from the illuminating stars than the PAHs, lend- 
ing support to the idea that PAHs are produced from 
the photo-evaporation of these very small grains. While 
more examples are needed for understanding the origin 
and fate of interstellar PAHs, it appears nonetheless that 
these macro molecules are released in the gas phase in the 
UV-illuniinated regions of the interstellar medium, i.e. in 
the diffuse clouds, in PDRs, etc. In those regions, the de- 
struction of the carbon skeleton is the main process lim- 
iting the smallest possible PAH size. It is likely that some 
carbon bearing molecules are released in the gas phase in 
the UV-illuminated regions, either as a secondary prod- 
uct of the evaporation of the dust particles giving rise to 
PAHs, or as products of the destruction of the PAH car- 
bon skeleton. 

5. Summary and Conclusions 

We have presented maps of the edge of the Horsehead neb- 
ula in rotational lines of excited H2, CO, C^^O and simple 
hydrocarbon molecules, CCH, C-C3H2 and C4H with 6" 
resolution. All the hydrocarbon maps are strikingly simi- 
lar to each other, and to the m id-IR emission mapped by 
ISOCAM ijAbergel et al.ll2003|) while we measured a 3 to 
5" offset between the hydrocarbon and H2 peaks. State-of- 
the-art chemical models fail to reproduce both the PDR 
hydrocarbon stratification and the absolute abundances of 
2 over 3 observed hydrocarbons. We have examined three 
hypotheses for improving the models, and we conclude 
that the most likely explanation is we are witnessing the 
fragmentation of PAHs in the intense far-UV radiation 
due to crOri. 

A detailed modeling of the chemistry including this 
new mechanism is beyond the scope of this paper. Indeed, 
such a modeling requires rates for both the growth (by 
addition of molecules or of carbon and hydrogen atoms) 
and the fragmentation of PAHs. This last item requires 
an accurate description of the fragmentation cascade of 



PAHs, in all their possible equilibrium states (ionized, neu- 
tral, partially or totally dehydrogenated, . . . ). Laboratory 
experiments such as t he ion cyclot ronic resonance cell 
PIRENEA in Toulouse l|.Toblinll20n3^ are key instruments 
in this perspective. In addition, the rate files used by the 
model need to be updated, with a special mention on the 
photo-dissociation rates of the simple carbon chains. A 
critical review of the role of neutral-neutral reactions in 
interstellar chemistry is also warranted. 
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